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ABSTRACT 

Pulsars in close, eccentric binary systems are usually assumed to have another neu- 
tron star as a companion. These double neutron star binaries have proven to be the 
best laboratories for experimental General Relativity and are the most secure candi- 
dates for gravitational wave interferometers. We present deep B, V, and R images of 
the field containing the eccentric binary pulsar system PSR B2303+46. We find a faint, 
blue object (B = 26.60 ± 0.09; (B - R) = -0.4 ± 0.2) coincident with the timing 
position. We suggest this object is the optical counterpart to the PSR B2303+46 sys- 
tem. The counterpart is too bright to reflect emission from the pulsar or a neutron star 
companion. Most likely, the companion of PSR B2303+46 is not a neutron star but a 
massive white dwarf. We show that the observations are consistent with a hot white 
dwarf companion (T e s ^5 x 10 4 K) with cooling age equal the characteristic age of the 
pulsar (tcooi — 30 Myr) and mass within the range set by timing observations and the 
Chandrasekhar mass (1.2 < Mq < 1.4 M & ). Given the eccentric orbit, the white dwarf 
must have formed before the neutron star, from what was originally the more massive 
star in the binary. Due to mass transfer, the originally less massive star could become 
sufficiently massive to end its life in a supernova explosion and form the radio pulsar. 
We constrain the mass of the pulsar to be in the range 1.24 < Mpsr < 1.44 M & . 

Subject headings: binaries: close — pulsars: individual (PSR B2303+46) — stars: evo- 
lution 
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1. Introduction 



2. PSR B2303+46 



Binaries in which a radio pulsar is in a close, ec- 
centric orbit around a compact companion have pro- 
vided the best astronomical tests of General Relativ- 
ity and the most accura te neu tron-star mass deter- 
minations (Taylor et al. 1992 ); are the most secure 
targets for gravitational wave interferometers such as 
LIGO (Abramovici et al. |1992| ); and may be the pro- 
genitors of the enigmatic gamma-ray bursts (Piran 
1997 ). Especially for the latter two topics, the na- 
ture of the compact companion is of vital importance. 
Fortunately, optical observations can distinguish the 
three possibilities: white dwarfs can be detected out 
to large distances, whereas neutron stars are rather 
dim and black holes not detectable at all. 

It has become general practice to identify as double 
neutron star system any high-eccentricity binary pul- 
sar for which the inferred companion mass is ~ 1.4 M Q 



(for a review, Van den Heuvel 1995). This is based on 
the evolutionary scenario for these systems. Briefly, 
it starts with an early-type binary. The primary (the 
more massive star) evolves first, transfers some mass 
to the secondary, explodes as a supernova, and forms 
a neutron star. Next, the secondary evolves and 
starts to transfer matter to the neutron star. The 
accretion causes "recycling" : spin-up and - in a way 
not well understood - a reduction in magnetic field 
strength. The mass transfer will be unstable, leading 
to a common-envelope phase. The orbit will necessar- 
ily be circularized; it can only become eccentric if the 
secondary is sufficiently massive to explode and form 
a neutron star in turn. (Otherwise, a white dwarf is 
left in a circular orbit, as observed for other binary 
pulsars.) Thus, if the observed radio pulsar can be 
shown to be the first-formed, recycled neutron star, 
the companion must be a neutron star as well. 

In four of the six presumed double-neutron star 
systems - PSR B1913+16, J1518+4904, B1534+12, 
and B2127+11C - the pulsars indeed appear recy- 
cled: the spin periods are shorter and inferred mag- 
netic fields weaker than for ordinary pulsars (tens of 
ms vs. ~ls, and ~ 10 10 G vs. ;> 10 12 G). However, 
the pulsars in the other two systems, PSR B1820— 11 
and B2303+46, show no clear sign of recycling. For 
the former, Phinney & Verbunt (1991) suggested that 
the companion was not a neutron star but a low-mass 
main-sequence star, and that the system will even- 
tually become a low-mass X-ray binary. The latter, 
PSR B2303+46, is the subject of this paper. 



PSR B2303+46 is in a 12.3d, highl y ecce ntric (e = 
0.66) orbit (Stokes, Taylor, & Dewey |l985|) . Perias- 
tron advance is observed, from which one infers a to- 
tal mass of the system M PSR + M c = 2.64 ± 0.05 M Q , 
as well as, in combination with the mass function, 
the limits M PS R < 1.44 Mq and M c > 1.2OM 



(Thorsett et al. |199J ; Thorsett & Chakrabarty |1998| ). 

The pulsar period (P — 1.06 s) and inferred mag- 
netic field strength (B ~ 8 x 10 11 G) allow the possi- 
bility of mild recycling. Given the small characteristic 
age, tpsr = P/2P ~ 30Myr, radio pulsations of the 
presumed neutron-star companion have been searched 
for, but with no success. Kulkarni (1988) carried out 
optical observations with the Palomar 200-inch, but 
found no counterpart down to R = 26 mag; this was 
seen as confirmation of the scenario outlined above. 

The pulsar parameters, however, are also consis- 
tent with those of ordinary pulsars. Thus, the pul- 
sar could have have formed after the companion com- 
pleted its evolution. If so, it may be the only neutron 
star in the binary, the companion being a white dwarf. 
This requires a twist to the evolutionary scenario, in 
which one starts with two stars with masses (slightly) 
below the critical mass, Af cr i t (~ 8Af Q ; Koester & 



Rcimers 1996), required to evolve to a neutron star. 
In due course, the primary evolves, transfers matter 
to the secondary, and forms a white dwarf. Now if 
the mass transfer increased the secondary mass be- 
yond Merit, h can explode and form a neutron star, 
resulting in a binary with an older white dwarf and a 
younger neutron star in a highly eccentric orbit. 

We were reminded of this possibility by Dr Wij- 
ers (1997, private communication), who wondered 
whether it could be verified observationally. At a dis- 
tance of 4.3 kpc (inferred from the pulsar dispersion 
measure) and for a cooling age of 30Myr (tpsr), a 
1.2 Mq white dwarf counterpart would have V ~ 25. 
This is excluded by the limit mentioned above, but, 
perhaps fortuitously, we had forgotten about this re- 
sult. Here, we report new, deeper optical observations 
with the Keck telescope, which show a possible coun- 
terpart to the PSR B2303+46 system. 

3. Optical Observations 

We imaged the field containing PSR B2303+46 
with the Low Resolution Imaging Spectrograph (Oke 
et al. 1995[ ) at the Keck II telescope, on the nights of 
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November 28 and 29, 1997 (UT). On the 28th, three 
600-s exposures were obtained in the R band, and two 
900-s exposures in B. On the 29th, one 450-s and four 
600-s in R, three 900-s in B, and five 600-s exposures 
in V were taken. All images were taken at airmass 
< 1.4. The skies were clear on the second night, but 
the first night was plagued by cirrus. 

The reduction was done as described by Kulkarni & 
van Kerkwijk ( fl99gj ) for the field of RX J0720.4-3125, 
which was observed on the same nights. For the pho- 
tometric calibration, we used Landolt fields: in B and 
R, the four listed in Table 1 of the above reference; 
in V, the first two only. We estimate uncertainties of 
<j0.02 mag in the zero points. 

For the astrometry, we selected from the USNO- 



A2.0 catalogue (Monet et al. |1998[ ) all 163 stars that 
overlapped with a 10-s R-band image. We measured 
their centroids, and corrected for instrumental distor- 
tion using a bi-cubic function determined by Cohen 
(private communication). With the plate scale known 
accurately, we fitted only for the zero-points in each 
coordinate and the position angle on the sky. Af- 
ter rejecting 7 outliers (residual larger than 0'.'8), the 
root-mean-square residuals were 0'.'20 in each coordi- 
nate. The astrometry was transferred to the stacked 
B, V, and R images using 28 transfer stars close to the 
pulsar position, solving again for rotation and zero- 
points. The rms residuals were <;0'.'04. 

Close to the timing position of PSR B2303+46, 
we found one faint, relatively blue object, hereafter 
star 1; see Figure [j] and Table [l| To verify whether the 
respective positions are consistent with each other, 
one has to take into account the measurement uncer- 
tainties (cti = 0'.'04, crpsR = 0'.'17 in each coordinate), 
as well as the extent to which the two positions are 
on the same astrometric system. The USNO-A2.0 
catalogue is tied to the International Celestial Ref- 
erence System as realized by the Tycho-based ACT 



catalogue (see Monet et al. |1998| ), while the DE200 
dynamical cphcmcris - on which the pulsar timing po- 
sition is based - is close to the ICRS as well (Folkncr 
et al. 1994). We expect the frame difference between 



DE200 and USNO-A2.0 to be considerably less than 
cpSR- Star 1 is offset by C'28 from the timing posi- 
tion, i.e., well within the 95% confidence radius of 0'.'4 
one infers from opsr alone. 

For the photometry, we used a simple point-spread 
function fitting method, which takes into account gra- 
dicnt s in the sky level (see Kulkarni & van Kerkwijk 
1998 ). First, we used aperture photometry on im- 



Table 1 
Photometry and Astrometry 



Star a.72000 

^ h m s 



<5.J2000 
) ( ° ' " ) 



B 

(mag) 



V 
(mag) 



R 

(mag) 



PSR 23 05 55.842 



23 05 55.869 
23 05 56.875 
23 05 59.665 
23 05 59.362 
23 05 57.242 
23 05 56.734 
23 05 56.574 
23 05 56.036 
23 05 56.682 



47 07 
47 07 
47 07 
47 07 
47 07 
47 07 
47 07 
47 07 
47 07 
47 08 



45.32 
45.30 
29.37 
43.54 
29.88 
51.89 
54.47 
35.96 
38.18 
03.55 



26.60( 9) 

19.84( 2) 

24.08( 2) 

22.79( 2) 

22.00( 2) 

25.18( 3) 

26.51( 8) 
26.73(10) 

26.43( 8) 



26.91(20) 
18.89( 2) 
22.64( 
21.16( 
20.62( 
23.46( 
24.74( 
26.39(11) 
26.20( 9) 



26.65(16) 

21.72( 2) 

20.10( 2) 

19.69( 2) 

22.21( 2) 

23.40( 2) 

25.99( 9) 

25.87( 8) 



Note. — The top line gives the 



timing position for 
PSR B2303+46 (Thorsett et al. |1993[ ); the uncertainties are 
0?017 and 0''17. Star 1 is the proposed optical counterpart; 
other stars are discussed in the text. The measurement un- 
certainties for the positions of optical objects are <0?004 and 
<0'/04. Possible uncertainties in the tie to astrometric systems 
are discussed in The uncertainties in the photometry are 
indicated by the numbers in brackets. Star 2 has no R-band 
magnitude, as it was overexposed in the R-band images. 



ages from the 29th to measure the instrumental mag- 
nitudes for relatively isolated, brighter "secondary" 
stars (stars 2-6 in B, V; 3-6 in R, because star 2 was 
overexposed; Fig. [l] and Table [I]). Next, we extracted 
from the stacked images 21 x 21 pixel (4'.'5 x 4'.'5) re- 
gions around stars 2-6, the candidate (star 1), and 
three other faint objects (stars 7-9). We fitted these 
to a 2-dimensional Gaussian on top of a plane with 
an arbitrary tilt, and determined the average FWHM 
for the secondary stars. We then refitted all objects 
keeping the FWHM fixed at the average, and used 
the amplitudes of the Gaussians to determine relative 
magnitudes. Finally, the difference with the aper- 
ture results for the secondary stars was used to calcu- 
late instrumental magnitudes for stars 1, 7, 8, and 9, 
and all magnitudes were calibrated using the solu- 
tion found from the Landolt stars. The results are 
listed in Table 13. In order to verify our procedures, 
we also determined B and R magnitudes for stars 1, 
7, 8, and 9 from the stacked images from the 28th and 
29th separately: these gave consistent results. 

4. A Massive White Dwarf Companion 

The probability that star 1 is a background ob- 
ject that happens to be within the 0.5arcsec 2 er- 
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Fig. 1. — Stacked B, V, and R images of the localization of PSR B2303+46. In each, the timing position is indicated 
with the tick marks (of size 1"5). Stars mentioned in the text are labeled below their image. 



ror region (95% confidence) of the timing position of 
PSR B2303+46 is about 5%, i.e., not particularly low. 
With B — R = -0.05 ± 0.18, however, star 1 is bluer 
than all other faint stars in the field. In our images, 
the bluest other objects have B — R~ 0.6 (e.g., stars 8 
and 9). Of these, there are only a few arcmin~ 2 and 
the chance coincidence probability is < 0.1%. The 
probability for an object as blue as star 1 is lower 
still, and, therefore, we believe star 1 is the optical 
counterpart to PSR B2303+46. 

It is unlikely that the optical emission is due to the 
pulsar or to a neutron-star companion. Thermal emis- 
sion from a neutron star could reproduce the colors, 
but it would be much too faint: the known sources 
have similar magnitudes, but are all nearby (for a 
recent compilation, see Mignani 1998| ). Nonthermal 
emission can lead to brighter sources, but only for 
young pulsars and generally with colors that are too 
red. This leads us to propose that the companion of 
PSR B2303+46 is a massive white dwarf, and that 
star 1 is its optical counterpart. 

In order to verify whether our observations are con- 
sistent with a white-dwarf companion, we need to es- 
timate the expected brightness. This is possible using 
cooling models for white dwarfs, provided we have es- 
timates for the white-dwarf mass, composition, and 
age, as well as for the distance and reddening. We 
will discuss these in turn. 

The mass of the white dwarf companion has a strict 
lower bound of 1.2 Mq (inferred from timing; §H). An 
equally strict upper bound of 1.4 Mq is set by the 
Chandrasekhar mass. 



The age of the white dwarf is the sum of £sn, the 
time that elapsed between the formation of the white 
dwarf and the supernova explosion, and £psr, the age 
of the pulsar (see §||). An upper limit to £sn is set by 
the total lifetime of an 8- M (M c r it) star, i.e., ^sn < 



40Myr (e.g., Schaller et al. 1996 ). An upper limit to 
ipSR is set by the characteristic age (but see §||), i.e., 
tpsR <; tpsr = 30Myr. 

The cooling age of the white dwarf may equal its 
actual age, or i coo i = <sn + £psr < 70 Myr. It is quite 
likely, however, that the white dwarf was reheated 
when, prior to the supernova explosion, the system 
went through a common-envelope phase (required to 
account for the current small orbital size). If so, igN 
is irrelevant, and i coo i — ipsR <J 30 Myr. 

The composition of massive white dwarfs is still un- 
certain, but observations of novae indicate that both 
C+O and O+Nc+Mg are possible^ (Starrfield |1989| ). 
Fortunately, the mass-radius relations are very simi- 
lar (as inferred from Hamada & Salpeter 1961), and 
hence so should the cooling tracks. For massive white 
dwarfs, effects related to the composition of the at- 
mosphere are expected to be small as well (Wood 
1995| ). For completeness, we note that a hydrogen 



atmosphere seems likely, as some hydrogen will have 
been accreted during the common-envelope phase. 

The distance towards PSR B2303+46 can be con- 
strained from the observed dispersion measure (DM) 
of 62cm~ 3 pc. Using the Galactic model for the dis- 



For some massive white dwar, 
composition (Provencal et al 
Fe white dwarfs M max 



asses and radii indicate Fe 
This is excluded h ere: fo r 
1.1 Mq (Hamada & Salpeter |l96l|). 
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One infers 
and Ar — 



tribution of free electrons of Taylor & Cordes ( 1993 ), 
we find that in the pulsar's direction (7 11 = 105?41, 
& 11 = —11? 93) the predicted DM is consistent with 
the observed one for any distance d > 2.5 kpc. There 
is no upper limit on the distance, as the predicted 
maximum DM (for objects well outside the electron 
layer) is 78 ± 20 cm -3 pc. At the lower limit, the pre- 
dicted DM is 41 ± 10cm _3 pc. 

The reddening along the line of sight, estimated 
from dust infrared emission, is Eb-v = 0.22 ± 0.03 
(Schlegel, Finkbeiner, & Davis [19981 ). 
A B = 0.95 ± 0.13, A v = 0.73 ± 0.10 
0.59 ± 0.08. For any likely identification, star 1 is 
well out of the Galactic plane, and the full reddening 
should be taken into account. 

With the above, we are in a position to estimate 
the cooling flux from a presumed white dwarf com- 
panion to PSR B2303+46. From the cooling tracks 
of Benvenuto & Althaus ( 1999| ), we find that af- 
ter 30Myr a 1.2 M Q C+O white dwarf has cooled 
down to Mboi = 6.5 and T cff = 5 x 10 4 K. (Simi- 
lar res ults are obtained extrapolating tracks of Wood 
| 1995| ). From the atmos pheric models of Bergeron, 
Wesemael, & Beauchamp ( 1995 ), we find for this tem- 
perature BCv — —4.5 and Rayleigh- Jeans-like colors 
(B - V) = -0.28, (V - R) = -0.14 (for a hy- 
drogen atmosphere; differences for other atmospheric 
compositions should be small). Thus, one expects 
M B = 10.7 and (B - R) = -0.42. 

For a cooling age of 70Myr, we find Mboi = 7.5, 
T off = 4 x 10 4 K, M B = 11.0, (B - R) = -0.42. For 
higher masses, predictions are harder to make, as no 
cooling tracks are available. Extrapolating, we expect 
more massive white dwarfs to have smaller radii, but 
be hotter at the same age. The likely net effect will 
be that they are brighter bolometrically, but fainter 
in the optical (showing the same color). A conserva- 
tive lower limit Mb < 13.4 is inferred by taking into 
account the change in radius only, from ~ 0.006 Rq 
for a 1.2 Mq white dwarf to ~ 0.002 Rq for one at the 



Chandrasekhar mass (Hamada & Salpeter 1961). 

In summary, we expect that a white-dwarf com- 
panion will have (B - R) = -0.42 and 13.4 > M B ^ 
10.7. It can be brighter only if i coo i is substantially 
smaller than rpg^, which we consider unlikely. The 
expected reddened color, B — R = — 0.06±0.05, is con- 
sistent with the observations. The expected distance 
modulus is 12.3 < B — Mb — Ab < 15, corresponding 
to a distance of 3 < d < 10 kpc, consistent with the 
lower limit set by the dispersion measure. 



5. Ramifications 

We presented Keck imaging of PSR B2303+46, in 
which we identified a faint blue object, star 1, coinci- 
dent with the precise timing position. We have shown 
that the companion to PSR B2303+46 could be a hot, 
massive white dwarf, with star 1 its optical counter- 
part. We have outlined how the evolution of a binary 
composed of two stars with masses close to but be- 
low M cr i t could lead to the formation of a system like 
PSR B2303+46. 

The system may help calibrate different chronome- 
ters. With a temperature measurement from UV ob- 
servations (HST time granted), cooling models will 
provide a lower limit to t coo \. This may allow a test 
of pulsar braking: if one finds t coo i 3> ^sn + Tpsr> 
the braking index n has to be the culprit. The brak- 
ing index enters via tpsr = P/(n — 1)P. Usually, 
n = 3 is assumed (as we have done above), which 
is valid for a dipole rotating in vacuo. However, this 
has not been verified observationally. Indeed, all mea- 
surements give n < 3, although these were for some 
young, fast pulsars, for which n < 3 is perhaps ex- 
pected (Melatos 1997). It would be interesting as well 
to find t coo \ <C tpsr: this would imply that pulsars 
do not have to be born with short spin periods. 

If the companion is indeed a white dwarf, both the 
lower and the upper limit to its mass are interesting. 
The upper limit, in combination with the total mass, 
corresponds to a lower limit to the mass of the pulsar. 
Combined with the upper limit set by timing (|2]), one 
finds 1.24 < A/psr < 1.44 M©. 

The lower limit of 1.2 Mq makes the companion in- 
teresting as a white dwarf. Masses > 1 M have also 
been inferred for about a dozen field white dwarfs, 
many discovered only recen tly from extre me UV sky 
surveys (Marsh et al |l997| ; Vennes et al. |l997| ; Fin- 
ley, Koester, & Basri 1997 ). It is not clear, however, 
whether these have been formed from massive stars, 
for the following reasons. First, quite a few, espe- 
cially the more massive ones, have strong magnetic 
fields, which is unusual; also, it makes the mass es- 
timates, which are based on line broadening, more 
uncertain (e.g., Schmidt et al. 1992 ; Ferrario, Vennes, 
& Wickramasinghe 1998). Second, most are hot and 
therefore young; combined with the short lifetimes of 
stars massive enough to form them, one would expect 
to find them in young star clusters, not in the field 
(Bergeron et al. 1991). Third, in at least one object, 
GD50, unexpected traces of helium have been found 
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in the spectrum, as well as evidence for a high rotation 
rate (Vennes, Bowyer, & Dupuis 199(f ). Fourth, the 
massive white dwarfs appear to form not just a tail of 
the distribution of white-dwarf masses, but rather a 



separate peak (Finley et al. 1997). These reasons have 
led to the speculation that at least some of these field 
massive white dwarfs are not the product of single- 
star evolution, but rather the result of mergers of 
two ordinary 0.6 M Q white dwarfs (references cited 
above). 

White dwarfs in young star clusters and binaries 
almost certainly are the product of massive stars, but 
these have masses up to ~ 1 M Q only (in NGC 2516 
[Koester & Reimers 1996[ | and Sirius [Gatewood & 
Gatewood |197§ ; Provencal et al. 1998]). In contrast, 
the massive white dwarf companion of PSR B2303+46 
has almost certainly descended from a massive star 
and it is undeniably massive - a statement that can 
be made given the exquisite precision of pulsar timing. 
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